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By Hsu Lo 

SUMMARY 


A short and simple method is presented for the determination of 
transient skin temperature of conical Bodies during short-time, high-speed 
flight. A differential equation is presented for this purpose, 
giving the fundamental relations Between the transient skin temperature 
and the flight history. For the heat— transfer coefficient and 
Boundary— layer temperature, which are needed in the differential 
equation, EBer*s experimental results for conical Bodies under super- 
sonic conditions are adapted and summarized in a convenient way. The 
method is applied first to flight at constant altitude to illustrate 
the effect of acceleration on transient skin temperature. The method is 
then applied to arbitrary flight. Several examples are given; for one 
example measured data are available and are in good agreement with the 
calculations . 


INTRODUCTION 


When air flows over a Body, the air immediately adjacent to the 
Body is Brought to rest By skin friction. As a result the air is 
heated to a higher temperature and, hence, heat exchange Between the 
air and the skin occurs. This phenomenon is generally termed 
"aerodynamic heating." 

At high speed the temperature increase of the air is very large 
and the aerodynamic-heating problem Becomes of great concern to 
designers. The problem is related to the characteristics of the 
Boundary layer and the local heat— transfer coefficient. In reference 1 
a method is given for the determination of the skin temperature at 
supersonic speed, which makes use of the formulas for heat— transfer 
coefficient and boundary— layer temperature derived for flat plates at 
subsonic speed. In reference 2 a different approach is made for the 
determination of skin temperature of a Body of revolution in supersonic 
flight. Both papers, however, deal with equilibrium skin temperature 
for steady-flight conditions at constant altitude. 
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In several German papers (reference 3, for instance) it is shown 
that for a short-time flight during which the speed and altitude vary 
with time, the transient skin temperature may be considerably lower 
than the equilibrium skin temperature. In the present paper, therefore, 
emphasis is given to the transient skin temperature, rather than the 
equilibrium skin temperature. A differential equation is presented for 
this purpose; and for the heat— transfer coefficient and boundary— layer 
temperature needed in the differential equation, Eber*s experimental 
results (reference 4) for conical bodies under supersonic conditions are 
adapted and summarized in a convenient form for immediate application. 

If, however, better experimental data become available, they can be 
adapted readily to the present method. 

In order to show the effect of acceleration on skin temperature the 
simpler problem of flight at constant altitude is treated first. More 
general flights are then discussed with several examples. In one example, 
data obtained from the Naval Research Laboratory, Washington, D. C., 
giving skin temperatures for the V— 2 missile 21 fired on March 7 , 19^-7 
at White Sands, N. Mex., in connection with upper atmosphere research 
(reference 5 ), are used for comparison with the calculated results. 


SYMBOIS 


a 

K 

c p 

c 

f 

g 

h 

h 1 


k 

2 

t 


acceleration, feet per second per second 
temperature recovery factor 

specific heat of air at constant pressure, Btu/(lb)(°F) 
specific heat of skin material, Btu/(lb)(°F) 
function of Tg and implicitly of t (see equation (l 6 )) 
acceleration due to gravity (taken as 32.2 ft/sec^) 
heat— transfer coefficient, Btu/(sec)(sq ft)(°F) 

reference heat— transfer coefficient corresponding 

to P = v (or 30 °), 2=1 foot, and H = Sea level, 

6 

Btu/(sec)(sq ft)(°F) 

thermal conductivity of air, Btu/(sec) (ft) (°F) 
characteristic length, feet 
time, seconds 


MCA TN No. 1725 


3 


r I 

G 

H 

J 

Ql 

0.2 

E 

Tl 

T a 

T e 

Ti 

T 

B 

T_ 


w 


Pi 


time when certain velocity is reached by uniformly 
accelerated body starting from rest, seconds 

correction factor for nose angle (see equation (l4)) 

correction factor for characteristic length (see 
equation (l4)) 

correction factor for altitude (see equation (l4) ) 
heat— absorption capacity of skin, Btu/(sq ft)(°F) 
altitude, feet 

mechanical equivalent of heat (taken as 778 ft— lb/Btu) 

heat flowing into skin due to skin friction, Btu/(sec) /(sq ft) 

heat lost through radiation, Btu/(sec)/(sq ft) 

Beynolds number 

free-stream temperature of air, °F absolute 
skin temperature at time t a , °F absolute 
equilibrium skin temperature, °F absolute 
initial temperature, °F absolute 
boundary— layer temperature, °F absolute 
skin temperature, °F absolute 
total or stagnation temperature, °F absolute 
free-stream velocity, feet per second 

specific weight of skin material, pounds per cubic foot 

total apex angle of conical nose, degrees or radians 

emissivity, assumed to be 0.4 in numerical examples 
of this paper (e = 1 for perfect black body) 

free-stream density of air, pounds per cubic foot 


MCA TN No. 1725 


free— stream density of air at sea level (taken as 
0.07657 Ib/cu ft) 


M- 


coefficient of viscosity of air, pounds per foot— second 


T 


skin thickness, feet 


ANALYSIS 


Fundamental Equations 


The skin temperature of a tody can he determined from a considera- 
tion of the amount of heat flowing into the skin and that lost ty the 
skin. The balance of the two is the heat absorbed or given up by the 
skin from which the change in skin temperature can be calculated. 

Heat flowing into skin ,— For heat flowing into the skin during a 
short-time, highr-speed flight, the most important factor to be cbn— 
sidered is the aerodynamic heating. Experimental results indicate 
that the heat flowing into the skin due to skin friction Qq (expressed 

in Btu/sec/sq ft) can be determined from the following empirical 
formula: 


where the boundary— layer temperature Tq and the heat— transfer 
coefficient h have been studied experimentally by Eber (reference 4) 
and are discussed subsequently. 

Heat received by the skin from other sources, like solar 
radiation, radiation from surrounding atmosphere, heat from interior 
or other parts of the body, and so forth, are not considered in this 
paper . 

Heat lost by the skin.— The heat lost by the skin can be heat 
lost through radiation, heat lost through artificial cooling, or 
heat lost to the interior or other parts of the body. The heat lost 
through radiation Q 2 (expressed in Btu/sec/sq ft) can be determined 
from the Stefan— Boltzmann formula: 



( 1 ) 



( 2 ) 


vhere the value of the emissivity e depends on the surface condition 
of the skin. For a perfect black body., 6=1. 
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The heat lost other than through radiation to the exterior is not 
considered in this paper. 

Heat balance e quation .— Equation (l) gives the rate of heat flow 
into the skin and equation (2) gives the rate of heat loss h y the skin. 
The difference of the two Qq — Qg is the heat left to heat the skin 
measured in Btu per square foot per second. During an interval dt, 
the total heat (Btu/sq ft) to he absorbed by the skin is therefore 

(Qq - q 2 ) at 


If the temperature rise of the skin during this interval dt is dTg, 
the heat absorbed by the skin can also be expressed by 


G dTg 


where G, the heat— absorption capacity of the skin (expressed in 
Btu/( sq ft) (°F) ), is the product of the specific heat of the skin 
material c, the specific weight of the skin material w, and the 
skin thickness t . From data in references 6 and 7* representative 
values of G for steel and alloys of aluminum and magnesium at 520° F 
are obtained and plotted against skin thickness in figure 1. At higher 
temperatures the values of G are increased at the rate of 0.256t’ 
per 100° F for steel and 0.0624t* per 100° F for aluminum where t* is 
the skin thickness in inches. For magnesium alloy the temperature effect 
can be neglected. 

Equating the preceding two expressions for the heat absorbed by 
the skin and substituting Qq and Q 2 from equations (l) and (2) 
gives 


G 



hTg 


+ 4.8 X 10 1 ^eT 0 ^ 


= hT 


B 


(3) 


which is the basic equation for transient skin temperature. 

T he simplified equation .- Equation (3) is a nonlinear differential 
equation with variable coefficients. A simplified equation can be 
obtained if the radiation loss can be neglected. The simplified 
equation is 


d.To 

+ hTg = hT B 


(4) 


G 


dt 
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As is shown subsequently for example 3 , the radiation loss ordinarily 
contributes only a small part to the transient skin temperature for 
short-time, high-speed flight. For such flights the use of equation (4) 
does not introduce any appreciable error; whereas a great deal of time 
and labor can be saved in the computation work. 


Equilibrium temperature .— If a body flies with a constant speed 
at a constant altitude for a sufficiently long time, the skin is heated 
to a temperature such that the heat flowing into the skin per second 
is for practical purposes the same as the heat lost, and the skin 
temperature becomes constant. The skin temperature then is called the 
equilibrium temperature T e , corresponding to that speed and altitude. 
The equilibrium temperature T e can be determined from equation (3) 


by dropping out the first term since 



Thus, 


4.8 X 10 _13 €T ^ + hT 0 = hT 


B 


(5) 


If the radiation loss can he neglected, the equilibrium skin tempera- 
ture is the same as the boundary— layer temperature. 

In equations (3), (4), and (5) the two parameters h and T3 

are needed before the equations can be solved. In the following section 
Eber f s experimental data on these two parameters are discussed and 
summarized. 


Eber f s Experimental Results 

In 1941 Eber (reference 4) made a series of wind-tunnel tests on 
cones of various vertex angles to determine the boundary— layer 
temperature T3 and heat— transfer coefficient h at high speeds. 

The results he obtained are adapted in this paper for the determination 
of skin temperature for conical bodies in flight. These results are 
discussed in the following sections and are summarized in a convenient 
form with certain simplifications. 

Boundary— layer temperature T3.— When the air stream is brought 
to rest isentropically, the temperature of the air is called stagnation 
temperature. When the air is brought to rest by skin friction, the 
process is usually not isentropic and the temperature of the air is 
lower than the stagnation temperature . In this case, if no heat flow 
across the skin occurs, the temperature of the air immediately adjacent 
to the skin is called the boundary— layer temperature T3 and is found 

to be closeiy related to the stagnation temperature. In the actual case 
for which heat exchange between the air and the skin occurs, the 
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boundary— layer temperature T B is only an artificial term used in the 

empirical equation (l) and may not necessarily he realized at any point 
in the actual boundary layer. 

The stagnation temperature Tq can be calculated theoretically 
from the following statement of Bernoulli* s equation: 


where and T-^ are the velocity and the temperature of the free 

air stream. Values of T-^ at various altitudes are given in 
reference 8 and are replotted in figure 2 of this paper. 

For constant value of c p , equation (6) becomes 


For the actual case, however, c p varies with temperature. In reference 1 
a correction for variable c p has been shown to lower the stagnation 
temperature considerably at high velocity (about 20 percent lower at 
Mach number 8). An exact solution of Bernoulli's equation for variable c p 

gives a set of stagnation-temperature curves for various altitudes and 
Mach numbers. In this paper, however, plotting the stagnation^ 
temperature rise T^ — Tq against free— stream velocity Vq is chosen 
instead of the stagnation temperature against Mach number. The result 
is a single curve for all altitudes. This curve, for which the 
derivation is given in appendix A, is plotted in figure 3, together 
with the stagnation-temperature rise for c p = 0.24 Btu per pound 

per degree Fahrenheit. For low speeds, no correction for variable c p 
is necessary. 

The boundary— layer temperature Tg is related to the stagnation 

temperature by an empirical factor K called the temperature recovery 
factor. The recovery factor E is defined as 



(6) 



( 8 ) 


Tt - Tl 


Eber’s experimental results show that E varies with the total vertex 
angle (B of the cone but is practically independent of velocity. The 
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variation of K with (3 is not very large. For (3 ranging from 20° 
to 50°, an average value K = O .89 can be used with a maximum error of 
about 2 percent. Therefore, 

T B ~ T 1 = °* 8 9 (T t - T x ) (9) 

Since the stagnation— temperature rise Trji — T^_ can be taken as a 
function of Yj_ only (fig. 3 ), the boundary-layer-temperature rise can 
also be taken as a function of V]_; therefore, there is a single curve 

for boundary— layer— temperature rise for all altitudes. This curve is 
shown in figure 4. For c p = 0.2k Btu per pound per degree Fahrenheit, 
the boundary— layer—- temperature rise can be given by the following 
expression: 


T B - T X = 



( 10 ) 


Heat-tranafer coefficient h .- The heat-transfer coefficient h 
can be determined from Eber’s empirical formula which is taken from 
reference 

h = ( 0.0071 + 0.015^ \/p)y E 0,8 (11) 

where the Reynolds number R, as defined by Eber, is equal to ppV^l/p.. 
Equation (ll) can be rewritten as 

h = (0.OO7I + 0.015 k\T P) p l°* 8 "^8 V l°* 8 ( 12 ) 


In the derivation of equation (ll)* Eher used the "boundary— layer 
temperature as the reference temperature for values of k and M** the 
thermal conductivity and coefficient of viscosity of air, respectively. 
For the characteristic length 2 , Eher used the entire length of the 
surface of the cone. Although different opinions exist regarding what 
is the correct length to he used as the characteristic lengthy when 
Eher 1 s experimental results are applied to actual flight conditions 
(for - instance, in reference 3 half the total length of the cone surface 
is used as the characteristic length), Eher l s original definition is 
followed in this paper. The value of h thus obtained represents an 
average value, instead of local value* of the heat— transfer coefficient. 
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A simplified, method for the evaluation of the factor k/p^*® is 
now given. For a given velocity and altitude, the boundary— layer 
temperature can be determined from figures 2 and 4, and values of k 
and p corresponding to this temperature can be obtained from 
reference 9 for temperatures below 2400° F absolute. The boundary— layer 
temperature was determined for various velocities at three different 
altitudes (sea level, 100,000 ft, and 190,000 ft) and the ratio k/p^*^ 
was calculated and plotted against velocity in figure 5. For boundary- 
layer temperature higher than 24-00° F absolute the curve is shown by 
dashes and is obtained by extrapolation. 

In figure 5 the two curves corresponding to 100,000 feet 
and 190,000 feet form two limiting values of k/pO-8. For any other 
altitudes from sea level to 370*000 feet, approximately, values 
of k/p0-8 against velocity fall within these two limiting curves. 

Since both limiting curves do not differ very much from the sea— level 
curve, using the sea— level curve for all altitudes up to 370,000 feet 
seems Justified. Therefore, k/p0«8 is a function of velocity only. 

In connection with this simplification, it should be kept in 
mind that values of k and p at high temperature are obtained by 
extrapolation and any effects of change of air composition at high 
altitudes on k and p are not considered. 

Equation (12) is now reduced to four factors which are respective 
functions of p, 2, altitude, and velocity. For convenience of 
computation, a reference heat— transfer coefficient h 1 , corresponding 

to 3=| (or 30°), 2=1 foot, and sea— level conditions, is computed 

from equation (12) and the results are plotted in figure 6. This reference 
heat— transfer coefficient h* is a function of velocity only. For 
other nose angles, characteristic lengths, and altitudes, the heat- 
transfer coefficient h is simply the reference heat— transfer coeffi- 
cient h* multiplied by the three correction factors Fp, Fj, 

and F h . Thus, 


h = F p FjF H h* 


(13) 


vhere 


0.0071 + 0.0154^/3 
0.0071 + 0.0154 \J1k/6 



F- 


H 


PA0 . 8 

PA 


( 14 ) 
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Values of Fp, Fj, and Fg are given in figures 1 , 8, and 9 , respectively. 
Values of Pi/p 0 j where p 0 is the free— stream air density at sea level, 
were taken from the NACA standard atmosphere table (reference 10) for 
altitudes below 65,000 feet and from tables V(a) and V(b) of reference 8 
for altitudes above 65,000 feet. 

Application of Eber 1 s results .— In the foregoing discussion, Eber*s 
experimental results on Tjj and h are represented by simple curves 
as functions of velocity and altitude. For any prescribed flight path 
for which the velocity and altitude are given as functions of time, values 
of Tj and h can be expressed as functions of time. Table 1 is 
prepared for this purpose and the operations in table 1 are self- 
explanatory. Now that Tjj and h are known as functions of time, 
equation (3) or (4) can be solved for the transient skin temperature and 
equation (5), for the equilibrium skin temperature. 

Certain facts should be borne in mind, however, in the application 
of Eber’ s work, particularly to high— alt i tude flight conditions. First, 
Eber’s experimental results were obtained over a limited range of 

Reynolds number of 2 x 10^ to 2 x 10^; for flight conditions where the 
actual Reynolds number falls outside of this range (which is usually 
the case for flight at high altitudes), extrapolation is necessary. 

Second, Eber’s experiments were carried out at low altitudes. At 
high altitudes where the air density is very low, the flight may 
enter into the slip— flow domain where the recovery factor, and 
possibly the heat— transfer coefficient, may be greatly Effected. 

(For illustrations, calculations were made to determine the Mach numbers 
and Reynolds numbers from the flight trajectory of the V— 2 missile used 
in example 1 of the paper. The results of these calculations are 
indicated by points in figure 10, where the curve dividing the two 
domains is taken from reference 11. For this particular example the 
flight enters the slip— flow domain as soon as an altitude of approxi- 
mately 130,000 feet is reached.) Third, in the evaluation of k and p. , 
the effects of change of air composition at high eltitudes are neglected 
and values of k and p. at temperature greater than 2h00° F absolute 
are obtained by extrapolation. Finally, the properties of the atmosphere 
at high altitudes are taken from the tentative tables of reference 8. 

With all these uncertainties the application of Eber’s results to high 
altitudes may introduce a large percentage error. However, since the 
heat— transfer coefficient at high altitude is small, the effect on skin 
temperature is not large. 

A general discussion of the validity of the extension of Eber’s 
experimental data to flight conditions is given in reference 12 together 
with an extensive study of the effect on skin temperature of various 
individual parameters, such as Mach number, altitude, emissivity, 
characteristic length, and conical nose angle, as well as the day and 
night conditions. 
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Solutions of Equations 


Solution of equation (3) .— Equation (3) is a nonlinear equation 
with variable coefficients, of the first order hut the fourth degree. 
If equation (3) is written as 


it is readily recognizable that equation (15) can be solved con- 
veniently by Runge and Kutta* s numerical method (reference 13) which 
is summarized in appendix B. In table 2 Runge and Kutta* s method is 
arranged in a suitable way for the solution of equation (15) . A 
convenient time interval At is first chosen. The smaller At is, 
the more accurate the results will be. The operation of this table 
starts on the first line and proceeds from left to right and then to 
the succeeding lines. The function f, corresponding to t Tg 

at its left, can be obtained from equation (l 6) . The final results 
of table 2 give the skin temperature at the end of the interval At. 
If the table is repeated, the skin temperature at the end of 2At, 
3At, and so forth, can be obtained. If the time interval chosen is 
not sm a ll, a correction can be made as explained in appendix B, 
where an illustrative example is also given. 

Solution of equation (4) .— Equation (k) is a linear differential 
equation of the first order provided the heat— absorption capacity G 
is considered to be independent of skin temperature . The general 
solution is 



where D, the constant of integration, cam be determined from the 
initial condition Tg = Ti at t = 0. The parameters h and Tp in 
equation (17) usually cannot be expressed in simple analytical terms. 
Equation (17)* therefore, has to be integrated numerically. Table 3 
is provided for this purpose. The time and labor required for carrying 
out the computations in table 3 are much less than those required 
for table 2. Sometimes the values in row 8 of table 3 become very large 
as the time increases. It is advantageous to stop at this point and start 
again from the beginning of table 3 with the skin temperature last obtained 
being used as the new initial temperature. 

Solution of equation ([?) .— Equation (5) is a simple quartic 
algebraic equation. By use of Eber’s experimental results on h 



(15) 


where 



(16) 


V 
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and T B , equation (5) can be solved for the equilihrium temperature T e 

for any altitude and velocity combination. For illustrative purposes, 

the variation of equilihrium temperature with velocity for several altitudes 

is shown in figure 11. The value of e used is 0.4. 


RESULTS AND DISCUSSION 


The method for the calculation of transient skin temperature, 
as discussed in the preceding sections, is applied first to flight at 
constant altitude to illustrate the influence of acceleration on skin 
temperature and then to arbitrary flight conditions. 


Flight at Constant Altitude 

For illustration of the effect of acceleration, the simplified 
equation (4) is used to determine the skin temperature for the following 
example . 

Assume that a body, starting from rest, is accelerated uniformly 
to 5000 feet per second. This speed is then maintained until the 
equilibrium temperature is reached. Subsequently, the body decelerates 
uniformly to zero velocity again. The question is, then, how does the skin 
temperature change with time during these three periods of flight. The 
body has a conical no3e angle of 30°, a conical surface length of 1 foot, 
and a skin heat— absorption capacity G = 0.6 Btu per square foot 
per degree Fahrenheit and is traveling at a constant altitude 
of 50,000 feet. 

Period of uniform acceleration .- The velocity of the body at any 
instant during this period is given by 


Y ± = at 


where a is the acceleration in feet per second per second. The initial 
condition is Tg = T]_ at t = 0. 

Since the velocity is given as a function of time, T^ and h 

can be determined from table 1 and the skin temperature from table 3* 
all as functions of time. The results are given in figure 12 for 
five different values of acceleration: a = 2g, 5Sj 10g, 50g, and °o. 

(For infinite acceleration the temperature— time curve is but a single 
point.) The dashed curve in figure 12 gives the skin temperature T a 

at the time t a when the body reaches 5000 feet per second. TJie 
larger the acceleration is, the lower the skin temperature will be. 
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In figure 13* the skin temperatures for a = 2g and lOg are 
plotted against together with the curve showing the variation 

of equilibrium temperature with velocity. Since radiation loss is 
neglected in this case, the equilibrium temperature is the same as the 
boundary— layer temperature. The temperature difference between the 
"transient* 1 and the "equilibrium" curves is the "temperature lag/* 
which is greater for larger accelerations, as shown in figure 13. 

Period of constant velocity .— If the body maintains its speed 
at 5000 feet per second after this velocity is reached, h/G and T B 

are no longer functions of time; and the solution of equation (17) becomes 

— tt 

T S = T B + De G (18) 


for t ^ t a where D, the constant of integration, is determined, from 
the condition Tg = T & when t = t a . The skin temperature at the time 
when the "body first reached the velocity of 5000 feet per second T a 
is different for different accelerations. 

Equation (l8) is solved for the five different accelerations 
a = 2g, 5g, lOg, 50gj and °°. Results are plotted in figure 14. 

(The temperature curves before the body reaches ^>000 ft/sec are taken 
from fig. 12.) The skin temperature approaches exponentially the 
equilibrium temperature T e , in this case equal to Tg corresponding 
to 5000 feet per second at the altitude of 50,000 feet. 

Period of uniform deceleration .— If after the body maintains 
5000 feet per second for a long time it starts to decelerate, the 
velocity at any instant is given by 


V]_ = 5000 + at 


where t starts when the body begins to decelerate. The initial 
condition is Tg = Tg at t = 0. Table 1 and table 3 again can be 
used to solve equation (17) and the results are plotted against velocity 
in figure 15 for two different decelerations a = -2g and — lOg. Again, 
the equilibrium-temperature curve is also given. The transient skin 
temperature in the case of deceleration is higher than the equilibrium 
temperature throughout the velocity range. The temperature lag is 
therefore on the adverse side. 


General Flight Conditions 

In the general case, the body is changing its altitude as well 
as its velocity. Three examples are given. In all cases, the 
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simplified equation (4) is used for the calculation of transient skin 
temperatures. In example (3) the more accurate method is also used and 
the relative importance of radiation loss is discussed. 

Example 1.— The V— 2 missile 21 fired on March 7* 1947 at 
White - Sands, N. Mex., had a conical vertex angle of 26°, a conical 
surface length of approximately 7 feet* and a skin of 0. 109— inch steel 
after a certain distance from the nose. The above specifications 
and the flight path of the V— 2 missile as shown in figure l6 were 
obtained from reference 5* The boundary— layer temperature and the 
hea-t>-transfer coefficient are determined by means of table 1 and are 
plotted against time in figure 17. The skin heat— absorption 
capacity G is obtained from figure 1 to be 0.47 6 Btu per square foot 
per degree Fahrenheit corresponding to 520° F absolute., which is 
conservative. The initial skin temperature is known to be 546° F absolute. 
Table 3 is used to determine the skin temperatures and the results are 
plotted against time in figure 18. 

In figure l8 the measured skin temperature for the steel skin 
of the same missile is also shown. The measured data were obtained 
from the Naval Research Laboratory. (See reference 5-) The agreement 
between the calculated and measured results is good. 

Example 2.— An arbitrary velocity— and— altitude diagram is assumed., 
as shown in figure 19* including descending path as well as ascending. 

The missile is assumed to have a nose angle of 30°* a characteristic 
length of b feet, and a skin heat— absorption capacity of 0.6 Btu per 
square foot per degree Fahrenheit. The boundary— layer temperature T3 

and heat— transfer coefficient h are determined by use of table 1 and 
are plotted against time in figure 20. Table 3 is then used to calculate 
the skin temperatures. The results are plotted in figure 21. The skin 
temperature during descent becomes higher and higher because of the 
greater density at the lower altitudes and the higher velocity as 
it comes down. 

Example 3.— In order to determine the relative importance of the 
radiation loss, the calculation of skin temperature for the missile 
in example 2 for the first 80 seconds is repeated, except the more 
accurate method is now used where the radiation loss is not neglected 
(the emissivity € is assumed to be 0.4) . Table 2 is used for this 
purpose and the results are plotted in figure 22 together with the 
results obtained from example 2 where the radiation loss is neglected. 

The discrepancy between the two is very s m a l l, approximately 3° F. 

In fact, for most missiles the radiation loss plays only a small 
part in the determination of skin temperature and the simplified method 
can be used to great advantage in saving time and labor. A simple 
criterion unfortunately does not exist for predicting when the 
radiation loss can be neglected. Generally speaking, if the skin 
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heat capacity is not too small, say not below 0.5 Btu per square foot 
per degree Fahrenheit, the skin temperature is not expected to be higher 
than 1000° F absolute and, if the radiation loss is not the dominating 
factor for a long period of time, the radiation loss can be neglected . 

If the final skin temperature is completely unknown, starting with the 
simplified method is advantageous . From the maximum skin temperature 
obtained, an estimate of the radiation loss and its effect on the skin 
temperature can be quickly made. For instance, the maximum skin temperature 
of the missile of example 2 is about 615° F absolute. This skin temperature 
is conservatively assumed to hold for the entire period of 80 seconds. The 
change of skin temperature due to radiation loss during this period is 
then 


AT S 


i(4.8 X l<r 13 6Tg|)t 

^(4.8 X 10 -13 X 0.4 X 615^) X 80 


= 3 . 7° F 


which is negligible. A similar calculation for example 1 indicates that 
for the first 65 seconds of flight the radiation loss affects the 
maximum skin temperature approximately 6° F. 

For flight conditions where the radiation loss is the dominating 
factor for a long period, as from the 80 th second to the 220th second 
in example 2, the radiation loss should be investigated. During this 
period, h is zero and equation ( 3 ) becomes 

dTo * _ 1 o 4 

G — £ + 4.8 x 10 -j - 3 €T q = 0 

dt s 


The general solution is 


Tg = ^14.4 x 10 -13 e£ + d) 1//3 (19) 

where D, the constant of integration, can be determined from the 
condition Tg = (Tg) when t = 80 and where fis) q q s ^ n 

temperature at the 80 th second. At the end of the 220th second the 
skin temperature can be calculated from equation ( 19 ) and is found to 
be lowered only 6° F. 
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CONCLUDING REMARKS 


1. A differential equation taking into account aerodynamic heating 
and body radiation is presented for the calculation of transient skin 
temperature for any prescribed flight history. Runge and Kutta's 
numerical method is recommended for the solution of the differential 
equation. 

2. A simplified differential equation which neglects body radiation 
is also given and can be used in many cases to great advantage in saving 
both time and labor. 

3. Eber’s experimental results on the boundary— layer temperature 
and heat— transfer coefficient, to be used in the differential equation, 
are summarized in a convenient way for immediate application. Tables 
and charts are also provided to facilitate the solution of the 
differential equation. 

A. The calculated skin temperature for a V— 2 missile is in good 
agreement with the measured data. 

5. The heat— absorption capacity of the skin has an important 
influence on transient skin temperature. The heat— absorption capacity 
is greater and, consequently, the temperature lag is larger if the 
skin is thicker, the material is denser, or the specific heat is higher. 

6. When the air is heating the skin, the temperature lag due to 
the heat capacity of the skin is in the favorable direction, that is, 
tends to lower the skin temperature. When the air is cooling the 
skin, the temperature lag is in the adverse direction, that is, tends 
to keep the skin at high temperature. 

7. Eber’s experimental work was conducted under certain limited 
conditions (short testing time, small temperature difference, limited 
Reynolds number, and so forth) . More refined experimental values for 
a wider range are desirable. 

8. Because the atmospheric properties at extremely high altitudes 
as given in NACA TN No. 1200 are tentative and, also, because at high 
altitudes the flight actually enters the "slip— flow" domain, a closer 
investigation of the problem at high altitudes is needed. 

9- For more accurate results, investigations of heat exchanges 
other than those considered in this paper are necessary . 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., December 8, 19^7 
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APPENDIX A 

STAGNATION-TEMPERATURE RISE FOR VARIABLE c p 


In the following discussion, a single curve for all altitudes is 
obtained for the calculation of stagnation- temperature rise at various 
velocities Vp when the specific heat c p is considered to he a 

function of temperature. The Bernoulli equation can he written as 


V dV + Jgc p dT = 0 


(Al) 


Integrate equation (Al) between state 1, the free-etream condition, 
and state 2 , where the velocity is zero. 


or 




(A2) 


Values of the integral in equation (A2) can he obtained from 
table 1 of reference 9 (see also page 58 , reference 9 ) for given values 
of Tp and Tp,, and the velocity Vp can be computed. A set of 

curves can thus be obtained. 

However, if TVp — Tp is plotted against Vp with Tp as a parameter, 

the set of curves practical ly fall into one single curve for all values 
of Tp ranging from 392° F absolute to 63 O 0 F absolute, corresponding 

to the minimum and maximum free— stream temperatures for altitudes from 
sea level to about 370*000 feet. In figure 23 two curves are shown, 
representing the two extreme conditions. Using a single curve for the 
calculation of stagnation -temperature rise is therefore justified. 

The boundary— layer— temperature rise Tp — Tp can be obtained by 

multiplying the stagnation— temperature rise by the recovery factor K. 

The result is a single curve similar to the curve of stagnation— - temperature 
rise but with all the ordinates decreased by the ratio K. Figure k 
shows the curve of boundary— layer— - temperature rise at various speeds Vp, 

corresponding to K = O. 89 . 
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APPENDIX B 

RTJNGE AND EUTTA’S METHOD OF NUMERICAL INTEGRATION 


Runge and Kutta’ s method of numerical integration can he applied 
to differential equations of the following type, provided that the 
initial condition is known. 


S = f(i,y> (Ei) 


The initial condition is that y = jq when x = xq. 

The derivation of Runge and Kutta 1 s method can be found in 
reference 13 . The following table is provided for the calculation. 


X 

y 

f 

f X Ax 

Operation 

*0 

xo 

f ( x o^yo) 


§(*L + 04 ) = • • • 

^-2 + ^3 = • • • 
Sum = . . . 

q — ^Sum = • • . 

X 

o 

^0 + 1*1 

f ( x 0 + lp y o + b-i) 

<1 2 

X n + — 
0 2 

^0 + 2^2 

f (*o + + ^ 2 ) 

*3 

Xq + Ax 

y 0 + 1 3 

f (x Q + Ax,y 0 + q 3 ) 

% 

Xp = Xq + Ax 

yi = yo + i 





For convenience of application to the present problem, Runge 
and Kutta* s method is arranged in the form of table 2 of this paper. 

A suitable interval At should be chosen. If the pair of initial 
values 1 q and T q are given, the values t* and T* at the end 

of the interval At can be obtained by carrying out the operations in 
table 2. The operations proceed from left to right of the first line 
and then the succeeding lines. The function f(t,Tg), corresponding 
to t and T to its left on the same line, can be obtained from 
equation (l6). If the pair of values t* and T* are used as 
initial values and the table repeated, another pair of values t n 
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and T !l can "be obtained, corresponding to the skin temperature at the 
end of 2At. Table 2 can he repeatedly used in this manner until the 
desired values of T are obtained. 

For illustration, table 2 is used to calculate the skin temperature 
of the missile in example 2. Assume that the skin temperature at the 
end of 20 th second is knovn to be 519 ° F absolute, and the skin 
temperature at the end of 24th second is to be determined. 

Choose At = 2 seconds. 


t 

T S 

h 

t b 

G 

f 

f X At 

Operation 

20 

519.0 

0.0103 

544 

0.6 

0.403 

q-^ = 0.806 

|(n + 14 ) = 1 - 6lS 

21 

519.4 

.0106 

568 

.6 

.834 

qg = 1.668 

Ip + 13 = 3.324 

21 

519.8 

.0106 

568 

.6 

.828 

q^ = I .656 

Siam = 4 . 942 

22 

520.7 

.0107 

590 

.6 

1.215 

q 4 = 2.430 

q = ASum = 1.644 

22 

520.6 

0.0107 

590 

0.6 

1.217 

q x = 2.434 

5(11 + 11 ^ - 3.299 

23 

521.8 

.0108 

619 

.6 

1.728 

q. 2 = 3.456 

Q.o + Q.o = 6.886 

23 

522.3 

.0108 

619 

.6 

1.715 

1 3 = 3.430 

Sum = 10.185 

24 

524.0 

.0109 

640 

.6 

2.082 

q^ = 4.164 

q = isum = 3.395 

24 

524.0 



I 1 




The skin temperature at the end of the 24th second is therefore 
524° F absolute. Notice that a constant value of G is used in the 
above computation. For a more accurate analysis the value of G 
should be based on Tg given. 

In order to improve the result, the following method, as given by 
Runge and Kutta, can be used for corrections. Instead of At = 2 seconds, 
table 3 is repeated with At = 4 seconds. 


t 

T S 

h 

t b 

G 

f 

f X At 

Operation 

20 

519.0 

0.0103 

544 

0.6 

0.403 

q 1 = 1.612 

§(<ll + qj = 4.976 

22 

519.8 

.0107 

590 

.6 

1.230 

q 2 = 4.920 

^2 + ^ = 9»70^ 

24 

521.6 

.0107 

590 

.6 

1.196 

q 3 = 4.784 

Sum = 14 . 680 

24 

523.8 

.0109 

64o 

.6 

2.085 

qj*. = 8.340 

on 

o^ 

CO 

11 

Jen 

II 

J 

24 

523.9 
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The correction is then given by 



523.9) = 0.007 


and. the corrected skin temperature at the end of the 24th second is 


T s = 524 + I = 524.007 


Bunge and Kutta 1 s method can also he applied to simultaneous 
equations of the type (Bl) or to differential equations of higher 
order. For details ^ see reference 13- 
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TABLE 1 


VARIATION OF Tg AND h WZTH TIME 
joiven P, 1, Fp (fig. 7), and F i (fig. 8)] 


© 

© 

G) 

© 

© 

© 

© 

© 

® 

t 

(sec) 

H 

(ft) 

(given) 

V 1 

(ft /sec) 
(given) 

T 1 

(°F ats.) 
(from 
fig. 2) 

Tb -Ti 
(°F) 

(from 

fig. 4) 

t b 

(°F ats.) 
© ♦ © 

h» 

(from 
fig. 6) 

f h 

(from 

fig. 9) 

Wh* 
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TABLE 2 


NUMERICAL SOLUTION OF MOKE ACCURATE TRANS IENT-SKIN— TEMPERATURE EQUATION 


^Initial conditions^ t q = . . • sec and Tq = • • . °F abs.; time interval At = • . . 


t 

T S 

h 

(from 
table l) 

T 

B 

(from 
table l) 

G 

f 

(equation 16 ) 

f X At 

Operation 

t 0 

T 0 





*1 

|(<ii + * 14 ) = * * * 

^2 + ^3 - • • • 

vf 

4l 

T n + — 
0 2 





*2 

vf 

m , q 2 

T ° + T 





*3 

Sum = . . . 
q = isum = . . . 

3 

tQ + At 

T 0 + ^3 



• 


<14 

a t’ = tQ + At 

a T« = T 0 + q. 







a If t 1 and T 1 are used as initial conditions and the above table repeated^ the skin temperature 
corresponding to t = 1 q + 2At can be obtained. " naca 


ro 

uo 
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TABLE 3 

NUMERICAL SOLUTION OF THE SIMPLIFIED 
TRANSIENT-SKIN— TEMPERATURE EQUATION 


[Given G and t£J 


Row 

Steps of solution 

n 

1 

2 

3 


• • • 

© 

t 






© 

h (from table l) 






(D 

h/G 






® 

|(@n + ©n-l) 

a 





© 

At = © n - © n _! 

a 





© 

X 

© 

a 





© 

1 

0 





© 

e® 






© 

T-g (from table l) 






© 

© X © X © 






© 

|(@ n + ©n-i) 

a 





© 

© x @ 

a 





© 

^Z© + Ti 
1 

T i 





© 

_ © 
s = w 







^hese spaces are to be left blank. 
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Figure 1.- Heat -absorption capacity of various metal skins at room temperature 
(520° F abs.). (Based on data from references 6 and 7.) 
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Free- stream temperature, 7/, °F abs . 


Figure 2 .- Free-stream temperature at various altitudes. (Data obtained from 

reference 8 .) 



Stagnation -temperature rise > CT r - T f ) 
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Figure 3.- Stagnation-temperature rise for both constant and variable Cp. 
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Figure 4.- The boundary -layer -temperature rise at various velocities 

(variable c p ). 
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Figure 5. - Variation of k/ n 


with velocity. 


; 

Reference heat-transfer coefficient } h , (sec)(sqft)(r) 
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0 2000 4000 6000 8000 10000 
Free-stream velocity , V t , ft /sec 

Figure 6.- Variation of reference heat -transfer coefficient with velocity, 
p = 30°; l - 1 foot; and H = Sea level. 
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Total vertex angle , deg 


Figure 7. - Correction factor for p . 
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Figure 8.- Correction factor for l. 
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Altitude, H , -Ft 


(a) From sea level to 160,000 feet. 



(b) From 160,000 feet to 320,000 feet. 
Figure 9.- Correction factor for altitude. 


J 




5 



Figure 10.- Different domains of flow conditions entered by V-2 missile 21, fired on March 7, 1947. 

(Curve dividing the two domains is taken from reference 11.) ^ 
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Ecjuili brium temperatu re , T e? °F abs. 



Figure 11.- Equilibrium temperature at various velocities and altitudes, e = 0.4. 



Figure 12.- Time history of skin temperature for uniformly accelerated body at constant altitude. 

H - 50,000 feet; p = 30°; 1 = 1 foot; G = 0.6 Btu/(sq ft)(°F); t a = Time to reach 5000 feet per second. 
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O (000 2000 3000 4000 5000 


Free- stream velocity } V t ? ft /sec 


Figure 13.- Comparison of transient and equilibrium temperature for 
uniformly accelerated body at altitude 50,000 feet, p = 30°; l = 1 foot; 
and G = 0.6 Btu/(sq ft)(°F). 
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Figure 14.- Time history of uniformly accelerated 
body at altitude 50,000 feet. Constant velocity after 
the body reaches 5000 feet per second, p = 30°; 
1=1 foot; and G = 0.6 Btu/(sq ft)(°F). 


Skin temperature ? T$ } °Fabs. 
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Figure 15.- Skin temperature of uniformly decelerated 
body at constant altitude of 50,000 feet, p = 30°; 
1=1 foot; and G = 0.6 Btu/(sq ft)(°F). 
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Figure 16.- Velocity-and-altitude diagram for V-2 
missile 21, fired on March 7, 1947. (Data obtained 
from reference 5.) 
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Figure 17.- Variation of Tg and h with time for 

example 1. 
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Figure 18.- Comparison of calculated and measured skin temperatures for V-2 
missile 21, fired on March 7, 1947. p = 26°; 1 = 1 feet; and t = 0.109 inch 
for steel. (Measured data obtained from reference 5.) 
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Figure 19.- Assumed velocity -and -altitude diagram for example 2. 
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Variation of Tg and h with time for example 2. 


Figure 20.- 


Skin temperature, ^ °Fabs. 



Figure 21.- Variation of skin temperature with time for example 2. p = 30°; l= 4 feet- and 

G = 0.6 Btu/(sq ft)(°F). 
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Figure 22.- Comparison of simplified method and the 
more accurate method in example 3. p = 30°; 
l = 4 feet; and G = 0.6 Btu/(sq ft)(°F). 


\ 


-p- 

ro 



\/, , ft /sec 


Figure 23.- Stagnation-temperature rise; variable Cp. 
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